STOCHASTIC 2-D NAVIER-STOKES EQUATION WITH 
ARTIFICIAL COMPRESSIBILITY 
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Abstract. In this paper we study the stochastic Navier-Stoltcs equation with 
artificial compressibiUty. The main results of this work are the existence and 
uniqueness theorem for strong solutions and the limit to incompressible flow. 
These results are obtained by utilizing a local monotonicity property of the 
sum of the Stokes operator and the nonlinearity. 



1. Introduction 

The stochastic Navier-Stokes equation is a weh accepted model for atmospheric, 
aero and ocean dynamics. Chandrasekhar [5] and Novikov [14] first studied the 
Navier-Stokes equation with external random forces. After that several approaches 
have been poposed, from the classic paper by Bensoussan and Temam [4] to some 
more recent results, e.g., by Bensoussan [3], by Flandoli and Gatarek [8] and by 
Menaldi and Sritharan [13]. 

This paper is concerned with the existence and uniqueness of strong solutions 
for the Stochastic 2-D Navier-Stokes equation with artificial compressibility in 
bounded domains. The concept of artificial compressibility was first introduced 
by Chorin [6, 7] and Temam [17, 18], in order to overcome the computational 
difficulties connected with the incompressibility constraint. Using the classical 
Sobolev compactness embedding and exploiting the classical Lions [12] method of 
fractional derivatives, Temam in his papers [17, 18] and in his book [19] (chapter 
3) proved the existence, uniqueness and convergence of the deterministic Navier- 
Stokes equation with artificial compressibility in bounded domains. 

In the rest of this section we formulate the abstract Navier-Stokes problem with 
artificial compressibility. We describe some standard well known results including 
the local monotonicity property of the Navier-Stokes operator. In Section 2 we 
establish certain new a priori estimates involving exponential weight for stochas- 
tic Navier-Stokes equation with artificial compressibility. These estimates play 
a fundamental role in the proof of existence and uniqueness of strong solutions 
proved in the second half of Section 2. The monotonicity argument used here is 
the generalization of the classical Minty-Browder method for dealing with local 
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monotonicity. This method was first used by Menaldi and Sritharan [13] and for 
multiphcative noise by Sritharan and Sundar [16]. For similar ideas see also Barbu 
and Sritharan [1, 2]. In the last part of Section 2 we discuss the convergence of the 
corresponding perturbed problem. Here the use of local monotonicity avoids the 
classical method based on compactness and thus the results apply to unbounded 
domains and hence the existence and the uniqueness as well as convergence to 
incompressible flow are new even in the deterministic case. 

2. Abstract Mathematical Framework and Local Monotonicity 

Let O C be a bounded domain (for the sake of simplicity) with smooth 
boundary, u the velocity and p the pressure fields. The Navier-Stokes problem 
(with Newtonian constitutive relationship and artificial compressible medium) can 

be written as follows 

9tU-j/Au+(u- V)u+i(Divu)u + Vp = f in L'^{0,T;M-\0)), 

2 (2.1) 

e9tP + Divu = in h'^{0,T-l?{O)), 

with the initial conditions 

u(0) = Uo in L2(0) and p(0) = po in h'^iO), (2.2) 

where e > is a vanishing parameter, uq belong to L^(0) = L^(C',]R^), v is the 
kinematic viscosity, p denotes pressure and is a scalar-valued function and the 
(force) field f is in L^(0, T; L^(C')). A solution (u,p) should belongs to the space 
L2(0, T; Hi(0) X L2(C')), with MKO) = H1(C', R^) and H-i(C) its dual space. The 
second equation in (2.1) is an artificial state equation of a slightly compressible 
medium and the extra term (l/2)(Divu)u is a stabilization term to handle the 
nonlinearity. The standard spaces used here are as follows: 
llJ(C') with the norm 

„ \ 1/2 

\Vvfdx) = ||v||, (2.3) 
o ' 

and L^(C') with the norm 

.1/2 

O 



vl^da;) = |v|. (2.4) 

Using the Gelfand triple Vl\{0) C l?-[0) C H-i(e') we may consider A or V 
as a linear map from 110(0) or \?{p) into the dual of 110(0) respectively. The 
inner product in the \? or \? is denoted by (•, •) and the induced duality by (•, •). 
Thus, for any u = (uj), v = (vj) and w = {wi) in Ho(0) and p in L^(0) we have 

(— z/Au, w) = 1/ >J / diUjdiWjdx, (2.5) 

(-Vp,w) = -V' / dipwidx= / pdiWidx = {p,Diyw) (2.6) 
i Jo Jo 

((u • V)v, w) = / UidiVjWjdx. (2.7) 
„ Jo 



and 
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It is clear that u i-^- Divu is a linear continuous operator from 110(0) into L^(C'). 
Next, an integration by parts and Holder inequality yields 

((u- V)v,w) = -((Divu)w,v) - ((u- V)w,v), (2.8) 
K(u-V)v,w)|<C5^||«,«;,||^,^„^,J|a,«,||^,^^,,,, (2.9) 



and in the right-hand-side wo can use L -norms to estimate the product UiV 



J- 



Lemma 2.1. For any real-valued smooth functions and ip with compact support 
in M?, the following hold: 

< lk5ivp||,J|V52VILi, (2.10) 

IK. <2||<^||^J|V<,. (2.11) 
Proof. The results stated above are classical and well known [10]. □ 
As in Temam [19] , chapter 3, the non- linear term is a trilinear continuous form 

on ml{0) X ej(e») x ej(0) 

S(u, V, w) := {B{u, v), w) := ^ ^ / [ui diVj wj - Ui diWj Vj]dx, (2.12) 
where 

B{u) = B{u, u) = [(u • V) i Div u]u. (2.13) 
We have the following lemmas. 

Lemma 2.2. Let u and w be in the spaces ]HIq(0,M^) and 'L'^{0,M?) respectively. 
Then the following estimate holds: 

|(S(u),w)| < 2||u||3/2 |u|i/2 ||w||l4(o,m=). (2.14) 
Proof. We observe that 

S(u, V, w) = / Ui{DiVj)wj dx + ^ / {div u)vjWj dx. 

Using the Holder inequality, 

2 

|6(U,V,W)| < ^ ||Mi||L4(ei,M2) |AVj|L2(ei,M2) ||Wj||L4(C>,K2) 

2 



+ M*'^ "k^(o,R^) ll^jlk''(ei,M2) ||wj||L4(e(,R2). 



2 

In Frobenius norm divergence can be estimated by gradient. Hence 

3 

|6(u,u,wr)| < -||u1|l4(o^r2) ||u|| ||w||l4(c,_r2). 

Using the equation (2.11) in Lemma 2.1 we get 

|(i3(u),w)| = |&(u,u,w)| < 2||uf /2 |u|i/2 ||w||l4(o,r2). 
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□ 

Lemma 2.3. Let u and v be in the space 110(0,11$.^). Then the following estimate 
holds: 

\{B(n) - B(v), u - v)| < ^ ||u - vf + |^ |u - vj^ Mt^^^,^ ■ (2.15) 

Proof. For any given u and v in 111^(0,8^) , wc have from (2.12) 

(B(u),u) = 0, (2.16) 

and 

(B(u,v),v) = ((u-V)v,v) + i((Divu)v,v) =0. (2.17) 
Then using (2.16) we obtain 

(B(u) -B(v),u- v) = (B(u),u) - (B(u),v) +(S(v),v) - (B(v),u) 

= -(B(u),v)-(S(v),u). (2.18) 

Now 

{B{n - v), v) = (((u - v) • V)(u - v) + i(Div(u - v))(u - v), v) 
= (B(u), v) + (B(v), v) - ((u • V)v + ^(Div u)v, v) 
-((vV)u+i(Divv)u,v). 
Hence using (2.16) and (2.17) we have 

(B(u - v), v) = (S(u), v) - ((v • V)u + ^ (Div v)u, v) 

= (B(u),v) + (B(v),u) - ((v V)u+ ^(Divv)u,v) 

- ((v V)v + i(Divv)v,u) 
= (B(u),v) + (B(v),u) 

- ((v- V)(u + v) + ^(Divv)(u + v),u + v) 

+ ((v- V)u+ ^(Divv)u,u) 

+ ((v • V)v + i(Divv)v,v). 

With the help of (2.16) and (2.17) last three terms of the right hand side van- 
ish.Thus 

(B(u - v), v) = (S(u), v) + (i3(v), u). (2.19) 
Thus (2.18) and (2.19) yield 

(B(u) - S(v), u - v) = -(^(u - v), v), (2.20) 
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which by Lemma 2.2 gives the estimate 

|(S(u) -S(v),u- v)| < 2||u- v||3/2 |u- v|i/2 ||v| 



where || • || and | • | denotes the norm in lll{0) and L^(C') respectively. 

Now using the fact that for any two real numbers a, b and any p, q > 1 with 

^ + ^ = 1, 

P Q ' 

ab<^ + — , 

p q 

we obtain the estimate (2.15). □ 

Notice that Hj(0) is continuously included in L''(0) = Ij'^{0,R^) and u i-^ 
(V • u)u is a (nonlinear) continuous mapping from ]HIq(C) into its dual ]HI~^(0). 
Hence the nonlinear operator B{-) can be considered as a map from M\{0) into 
the space H^^(C') fl L'''/^(C', M^). Then combination of previous lemmas yield the 
following local monotonicity property. 

Lemma 2.4. For a given r > 0, let us denote by the (closed) U^-ball in 

= {v e H1(0,M2) . ||v||^,^^ ^,, < r}, (2.21) 

then the nonlinear operator u i— > Aa + B{\x) := —vAn + [(u • V) + (1/2) Divu]u 
is monotone in the convex ball B^. i.e., 

(Aw, w) + (B(u) - i3(v),w) + ^ |wp > - ||wf , (2.22) 
Vu e Hj(e',IR2), V e Br and w = u - V. 
Proof. First, it is clear that 

(Aw, w) = i/||w||^, 

and the equation(2.15) yields 

(S(u)-i?(v),w)>--||w|p-^|wp. 
Summing these equations up we get the desired result (2.22). □ 



3. Stochastic 2-D Navier-Stokes Equation with 
Artificial Compressibility 

Let us consider the Navier-Stokes equation subject to a random (Gaussian) term 
i.e., the forcing field f has a mean value still denoted by f and a noise denoted by 
G. We can write (to simplify notation we use time-invariant forces) f{t) = i{x,t) 
and the noise process G{t) = G{x,t) as a series dGk = J2kSk{x,t)dwk{t), where 
g = (gi,g2,''') and w = {^vi,W2, ■ ■ ■) are regarded as i'^-valued functions in 
X and t respectively. The stochastic noise process represented by g{t)d'w{t) = 
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J2k Sk{x,t)dwk{t,co) is normal distributed in HI with a trace-class co- variance op- 
erator denoted by = (t) and given by 

' (g'(t)u,v) = J2iSk{t),u) (gfc(t),v), 

(3.1) 

Tr(g^(t))=^|gfe(t)P<oo. 

k 

We interpret the stochastic Navier-Stokes equations as an Ito stochastic equa- 
tions in variational form 

' d(u(i), v) + {-iyAu{t) + [{u{t) • V) -h i Div u{t)]u{t) 

+ Vp(i),v)di = (f,v)di + ^(gfe,v)d«;fc(t), (3.2) 

fe 

. (ep(t)+Divu(i),(7) =0, 
in (0, T), with the initial condition 

(u(0),v) = (uo,v) and {p{0),q) = {po,q), (3.3) 

for any v in the space Hj(C') and any q in L^(0). 

A finite-dimensional (Galerkin) approximation of the stochastic Navier-Stokes 
equation can be defined as follows. Let {ei,e2, . . .} be a complete orthonormal 
system (i.e., a basis) in the Hilbert space L^(0) belonging to the space M.q{0) 
(and L4(e))). Denote by (£)) the n-dimensional subspace of ]L^{0) and Hj(C) 
of all linear combinations of the first n elements {ei, 62, . . . ,e„}. Also denote by 
L2 (O) V • hliO) the image of V. 

Consider the following stochastic ODE in R" 

' d(u"(t),v) + i-uAu^it) + [(u"(i) • V) + iDivu"(t)]u"(0 

+ Vp"(0,v)di = (f,v)di + ^(gfe,v)dwfe(t), (3.4) 

fe 

^ {ep"{t) + BiYu" it), q) = 0, 
in (0,r), with the initial condition 

(u(0),v) = (uo,v), (3.5) 

for any v in the space IL^(O) and q in L^(0). The coefficients involved are locally 

Lipschitz and wc need some a priori estimate to show global existence of a solution 
u"{t) as an adapted process in the space C'^(0, T, L^j(C')). 

Proposition 3.1 (energy estimate). Under the above mathemafAcal setting let 

f e l2(0,T;L2(O)), g e L\0,T; £2(1.^0))), uq G l2(0), po e l2(0). (3.6) 

Let u"(t) be an adapted process in C°(0,T, EI„) which solves the stochastic ODE 
(3.4)- Then we have the energy equality 

' d[|u"(i)|2 + £|p"(t)|2] + 1v |Vu"(t)|2dt 

= [2(f(t),u"(t))+Tr(g2(t))]di + 2^(gfe(i),u"(i))d«;fe(i), (3-^) 
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which yields the following estimate for any S > 

cT 



£;{|u"(t)|2+e|p"(f)|2}e-^* + 2!. / £;{|Vu"(i)|2}e-^*dt 

Jo 



< 



|u(0)p + £b(0)p + 1 [i|f(t)|2 + Tr(g2(0)]e-^*dt, 



(3i 



0<t<T 



for any < t < T. Moreover, if we suppose 

f e LP(0,T;L2(0)), g G LP(0,T;£2(L2(0))) (3.9) 

then we also have 

' E{jupJ\u-{t)\P +e\p^{t)\P]e-'' 

+ VV f |Vu"(i)nu"(i)|P-2e-'5*di} < |u(0)|P (3.10) 

Jo 

+ £|p(0)|P + C5,p,T / [|f(t)|P + Tr(g2(t))P/2]e-^*dt, 

JO 

/or some constant Cs,p,T depending only on 6 > 0, e > 0, l<p<oo and T > 0. 
Proof. Prom (3.4) we notice that, 

rf(u"(i),u"(i)) + (-z/Au"(i),u"(t))dt + ((u"(i) • V)u"(t),u"(i))di 

+ ^(u"(t) Divu"(t),u"(t))dt + (Vp"(t),u"(t))dt 

= (f(t),u"(0)di + ^(gfe(t),u"(t))dii;fc(t). (3.11) 

k 

It is clear that 

(-j.Au"(t),u"(i)) = i/|Vu"(t)|2, 
and the equation (2.8) yields 

((u"(i) . V)u"(i),u"(i)) = -(u"(i)Divu"(i),u"(<)) - ((u"(t) • V)u"(t),u"(t)). 
Hence 

((u"(t) • V)u"(t),u"(t)) + i(u"(t)Divu"(t),u"(t)) = 0. (3.12) 
Using the equation ep"{t) + Divu"(t) = we get from (2.6) 

(Vp"(t),u"(t)) = -(p"(i),Divu"(t)) = -(p"(t),-£p"(f)) = I ^b"(i)p. 

(3.13) 

Combining all the above results one can get from (3.11) 

^ d|u"(t)|2 + Z.|Vu"(t)|2dt + I d|p"(t)|2 



(f (i), u"(i))di + I Tr{g\t))dt + J2iSk{t), u^{t))dwk{t). 



k 



Rearranging the terms we get the desired energy equality (3.7). 
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Next, we calculate the stochastic differential of the process 
F{t) := [\u-{tr + s\p-{t)\']e-'' 

to get 

dF(t) =e-^*d[|u"(t)p+e|p"(t)|2] -5[|u"(t)|2+e|p»(t)|2]e-^*dt. 
Using the energy equality (3.7) we have 

dF(i) = - 2ry| Vu"(t)|2e-''*di + 2 (f (t), u"(t))e-''*di + Tr{g'^{t))e-^*dt 

+ 2Y,{Sk{t),n"{t))e-''dwkit) - d[\u^it)f + e\p^{t)f]e-''dt. (3.14) 

fc 

Now using the inequality 

2ab < 6a^ + ^b^ 





on 2(f(t),u"(t)) we have 

V /I - LI V ;i ' 1^ V ;i J ' ^ 



2(f(t),u"(t)) < 5|u"(i)P + ^|f(t)|2 < <5[|u"(i)P +eb"(t)|2] + -\m\\ 



Then (3.14) yields 

dF{t) = - 2zy|Vu"(t)|2e-**dt + (5[|u"(i)|2 + e\p'\t)\^]e'^*dt 

+ ^m\^e-''dt + Tr(g2(t))e-^*dt + 2 ^(gfe(t), u"(i))e-^*d«;fc(t) 

k 

-<5[|u"(i)|2+£|p"(t)|2]e-**di. 
Rearranging the terms we have 

d[{|u"(i)|2 + £|p"(t)|2}e-'5*] + 2i/|Vu"(t)|2e-**dt 



< 



hm\^e-''dt + Tr(g2(t))e-^*di + 2 V(g;t(i), u"(i))e-^*dw;fc(t). 
Next we integrate in [0, T] to get 

[|u"(i)P +e|p"(t)|']e-^* + 21/ r |Vu"(t)pe-**dt 

JO 

< |u(0)|^+£|p(0)|^ + ^^ [^mf+Trig\t))]e-''dt 

+ 2V r {gk{t),vi^{t))e-''dwk{t). 
k -^0 

Finally taking mathematical expectation and keeping in mind the expectation of 
a stochastic integral is zero, we have the desired result (3.8). 
Similarly, consider 

G{t) := [|u"(t)|P + £|p"(t)|P]e-^* 
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and use Ito calculus. Here we check that its stochastic differential satisfies 
dG{t) = - .5[|u"(t)|P + e\p"it)\P]e-^*dt 

+ ||u"(t)|P-2{d[|u"(t)|2+e|p"(t)|2]}e-^* 

Using the energy equality (3.7) we get 

dG{t) + J[|u"(t)|P + £|p"(t)|P]e-''*dt 

= ||u"(t)r'{-2H|u"(t)|rdt+[2(f(t),u"(t))+Tr(g2(0)]df 

k 

+ P(P_i)|u"(t)|P-4[4^(g,(t),u"Wfdt]e-^*. 



k 

Simplification and rearrangement of the terms in the above equation yields 
dG{t) + z/p||u"(i)|Hu"(t)|P-2e-''*dt + <5[|u"(i)|P + £|p"(i)|P]e-**dt 

= |u"(t)|P-2 [p m, u"(t)) + ^ Tviem e-'^dt 

+ vY.{Ek{t),u-m^''it)r^e-''dwk{t). (3.15) 

k 

Now in the first and the second terms on the right hand side of the above equation 
we apply the following elementary inequality 

Aab<l^ + «, 

p q 

where 1 + 1 = 1, A = a/3 > 0, ab > 0. 
Then choosing 



A = p, a = 1 , /3=(— . 

we get 



p (f (i), u"(i))|u"(t)|P-2 < aPM^ + K^!M^ 

= C,,p |f(f)|P + ^|u"(t)|P, 

where the constant C^^p > depends only on 5 > and 1 < p < oo. 
Similarly with proper choices of a and /?, we can prove that 

^TV(g2(i))|u"(0|P-^ < C,,p Tr(g2(i))« + ^|u"(t)|P. 
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Then (3.15) yields 

dG{t) + iyp\\u'\t)f\u^{t)\P-'^e-^'dt 

<C^,p[|f(i)r+Tr(g2(t))i]e-^*dt 

+pY.(^k{t),u-mn''{t)r^^~''dwk{t). (3.16) 

fe 

Integrating the stochastic differential (3.16),then taking the sup norm in [0, T] and 
finally taking the mathematical expectation we have 

e\ sup [|u"(t)|P + £b"(t)|P]e-** 

^ 0<t<T 

+ pz/^^|Vu"(t)nu"(t)|P-"e-^*di} < |u(0)|P 

+ £|p(0)|P + Q,p,T r [|f(i)r + Tr(g2(i))P/2]e-^*dt 
Jo 



+ 



pE{ sup /V(gfc(s),u"(s))|u"(s)|P-2e-*M«;fe(s)||. (3.17) 

0<t<T Jo , I-" 



fe 

By means of martingale inequality, we deduce 



e{ sup /*^(gfc(s),u"(s))|u"(s)|P-2e-^M^i;fc(s)|| 

0<t<T Jo ^ I 

<Ce{[ [(gfc(t),u"(i))|u"(t)|P-2e-^f di)'^'} 

<C e{(^£ Tr(g2(f))|u"(f)|2p-2e-2**dt) '^'} 

<C7i;( sup (|u"(t)|P-ie-**/p')f / TV(g2(i))e-2^*/Pdt)'^'| 

0<t<T ^Jo ^ J 

<^E{ sup (|u"(t)|Pe-**)} + C5,p,TM rTr(g2(t))P/2e-^*di}, 

^ 0<t<T Jo 

(3.18) 

where the constant Cs,p,T depends only on(5>0, l<p<oo and T > 0. 

Using (3.18) in (3.17) we get the desired estimate (3.10). □ 

Now we deal with the existence and uniqueness of the SPDE and its finite- 
dimensional approximation. 

Proposition 3.2 (uniqueness). Let u be a solution of the stochastic Navier-Stokes 

equation (SPDE) (3.2) with the regularity 

(ue L\n; C0(0, T; h^{0)) D L^{0, T; MI(0))), 

\ueh'^{nxOx{o,T)), L2(n X o X (o,r)), 



STOCHASTIC NAVIER-STOKES EQUATION WITH ARTIFICIAL COMPRESSIBILITY 11 



and let the data f , g, uq and po satisfy the condition 

( f e L2(0,T;H-1(0)), g e L2(0,T;4(L2(0))), 



(3.20) 



7/v in L^ifl; C°(0, T, L2(e'))nL2(0, T, Hj(0))) is anoi/ier solution of the stochastic 
Navier-Stokes equation (3.2), then 

rt 

' < 



27 
7^ 



4 

L^(0) 



|u(s)r. ds 



(3.21) 



[|u(<)-v(i)P+£b(i)-?(i)l'] exp^ 

< |u(0)-v(0)|2+£b(0)-g(0)|2, 

probability 1 /or any < t < T and e > 0. 

Proof. Indeed if u and v are two solutions then w = v — u solves the deterministic 
equation 

atw(t) - vA^{t) + V{q{t) - p{t)) = B{u) - B{v) in L'^{0, T; M-\0)), 

with B{u) = [(u- V) + i Div u]u. Notice that actually p and q are better processes, 
they belong to the space L2(0; (7^(0, T; L2(0))). 
Next, setting 

rit) :='^j\HsX,^^ds 

we have 

d(w(t),2e-''Ww(t)) - !/(Aw(t),2e-''(*)w(t))di 

+ {W{q{t)-p{t)),2e-^^'^w{t))dt 
= {B{u) - B(v),2e-''(*W(t))dt. 

Using (2.6) we get 

e-''Wd(|w(i)|2)+2i/e-'^(*)||w(t)||2dt + 2e-'^(*)(p(t)-g(t),Divw(t))dt 

= 2e-'-(*)(B(u) - B(v),w(i))di. (3.22) 
Since edt{q{t) - p{t)) + Div w(t) = 0, Lemma 2.3 and (3.22) yield 



\y,{t)f+e\q{t)-p{tt 

-2!ye-''W||w(i)||2dt - 2e-''W(B(v) - B(u), w(t))dt 

27 



< -2;.e-'-W||w(i)||'dt + 2e-'-(*)[^||w(t)||2 + _i|w(t)n|u(t)||^^,^^^ 
= -i/e-'-W ||w(t)||2dt + r(t)e-''W|w(t)|2dt 

< -i/e-'-W ||w(t)f di + r(i)e-'^(*) [|w(<)|2 + e|g(f) - p{t)\A dt. 



Hence 



-T{t) 



{Mt)?+e\q{t)-p{t)\''] 



< 0. 



Hence, integrating in t, we deduce (3.21), with probability 1. 



□ 
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Each solution u in the space L2(f]; L°°(0, T; 0-1(0)) n L2(0, T; Hj(C'))) of the 
stochastic Navicr-Stokes equation actually belongs to a better space, namely the 
space L2(n;C°(0,T;L2(C>)) nL4(0 x (0,T))) in 2-D, O C M^. Thus in 2-D, the 
uniqueness holds in the space L'^{n;L'^{0,T;B.-'^{O))). 

If a given adapted process u in L2(f2; L°°(0, T; h'^{O))nL'^{0, T; Hj(C»))) satisfies 

d(u(i), v) = (F(i), ^r)dt + {git), w)dw{t), (3.23) 

for any function v in M.q{0) and some functions F in L^{0,T;M~^{O)) and g in 
L^(0, T; £2(L^(C'))), then we can find a version of u (which is still denoted by u) 
in L2(f];C°(0,T;L2(e)))) satisfying the energy equality 

d|u(t)|2 = [2(F(t),u(t)) +Tr(g2(<)]di + 2(g(t),u(t))d«;(i) (3.24) 

see e.g. Gyongy and Krylov [9], Pardoux [15]. 

Definition 3.3. {Strong Solution) A strong solution u is defined on a given 
probability space {n, E, Et, M) as a L2(f2;L°°(0,T;L2(C))) n L2(0,r;Hi(O)) n 
C"(0, T; L^(C'))) valued function which satisfies the stochastic Navier-Stokes equa- 
tion (3.2) in the weak sense and also the energy inequality 

E{ sup [\u{t)\P + e\p{t)\P]e-^'+piy r \Vu{t)f\u{t)\^-^e-^*dt} 

0<t<T Jo 

< |u(0)|P + £b(0)|P + C5,p,T r [mi" + Tr(g'(i))P/'] e-**di, 

Jo 

where the constant Cs,p,T depends only on(5>0, e>0, l<p<oo and T > 0. 
Proposition 3.4 (2-D existence). Let f, g and Uq be such that 

r f e LP(o, T; m-\o)), g e lp(o, T; 4(l2(0))), 
\ uo e l2(o), po e l2(o). 

/or some p > 4. Then there is adapted processes u(i, x, w) ami p{t, x, u) with the 
regularity 

rueLP(fi;C0(0,T;L2(O)))nL2(O;L2(0,T;Hj(O))), 
\p, peL2(O;L2(0,T;L2(O))) 

satisfying the stochastic Navier-Stokes equation (3.2) and the a priori bound (3.10) 
/or every e > 0. 

Proof. Denoting 

F{u) := Au + B{u) - f := -uAu + [(u • V) -|- (1/2) Divuju - f 

we have 

du"(t) + F(u"(i))di + Vp"(t)dt = g{t)dw{t). 

Then using the a priori estimate (3.10), it follows from the Banach-Alaoglu theorem 
that along a subsequence, the Galerkin approximations {u"} have the following 
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limits: 



u" — >u weakly star in V'{Q:L°°{0,T;h'^{O))) nL^{n;L'^{0,T;W.l{O))), 

p" — >p weakly in L^{n;L^(0,T;h^(O))), 

i^(u") — >Fo weakly in L2(f];L2(0,r;H-i(C'))), 

where u has the Ito differential 

du{t) + Fo{t)dt + yp{t)dt = g{t)dw{t) in {n;L'^{0,T;m-\O))) 

and the energy equality holds, i.e., 

d[|u(f)|2 + s\p{t)\^] + 2{Fo{t),u{t))dt = Tr(g2(t))dt + 2{g{t),xi{t))dw{t). 

Now, for any adapted process v{t,x,u}) in L°°((0,T) x f2;L^(0)), we define 



27 /■* 

r{t,(j) := -3 / ||v(s,.,c 
Jo 



\\\ ds 



as an adapted, continuous (and bounded in u)) real- valued process in [0,T]. Then 
from the energy equality 

d[e-''W{|u"(t)|2+e|p"(t)|2}] +e-''W(2F(u"(t))+r(t)u"(i),u"(t))dt 
+ £|p"(t)|2r(t)e-''Wdi 

= TV(g2(t))e-''Wdi + 2(g(i), u"(i))e-''(*Mw;(i). 

Integrating between < t <T and taking the mathematical expectation we have 



E 



e-'-(^){|u"(r)|2+e|p«(T)|2}-|u 



"^'^^|2-e|p"(0)|2 



+ E 



■E 



f e-''W(2i?^(u"(i)) + r(Ou"(t), u"(t))dt 
Jo 



Jo 



|p"(t)pf(t)e-''(*Mi 



= E 



I TV(g2(t))e 
Jo 



Considering the fact that the initial conditions u"(0) and p"(0) converge to u(0) 
and p(0) respectively in L^, and the lower-semi-continuity of the L^-norm, we 
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deduce 



lim inf E 



-f 

Jo 



e-'-(*)(2i^(u"(t)) + r(t)u"(t), u^{t)}dt 



= lim inf E 



> E 



E 



e-''(^){|u"(T)|2 + £|p"(T)|2} - |u"(0)p - eb"(0)|2 

+ e [ |p"(t)|V(Oe-"(*Mt- / Tr(g2(f))e-'^(*Mt 
Jo Jo 

e-(^){|u(T)p + e\pm - |u(0)p - e\pm' 

r |p(i)|2f(t)e-'-(*Mi - r Tr(g2(f))e-'^(*M< 
Jo Jo 



+ e 



[ e-'-W(2J^o(t) + r(t)u(0, u{t))dt 
Jo 



Next, by monotonicity on L*-balls, i.e. by Lemma 2.4, we have 



2E 



f 

Jo 



e-'^(*)(i^(u"(i)) - i^(v(i)),u"(t) - v(i))dt 



+ E 



e-''«r(t)|u"(t)- v(t)|2di 



> 0. 



Rearranging the terms we find 



E 



Jo 



> E I e-''W(2F(u"(i)) +r(t)u"(i),v(t) - u"(i))dt 





Taking limit in n , we get 



E 



J 

Jo 



> E 



e-'-W(2F(v(t)) +r(Ov(t),v(t) - u(t))dt 



e-'-(*)(2Fo(i) +r(t)u(t),v(t) - u(i))di 



Now wc take v := u + Aw with A > and w is an adapted process in 
L4(f7; L°°(0, T; L2(C'))) n T; HjCO))). 

Then we have 

XE 



f e-'-W(2F(u(i) + Aw(i)) - 2Fo(t), w(t))dt 
Jo 



+ X^E 



J 

Jo 



e-''(*V(t)|w(i)|2di 



> 0. 



Dividing by A on both sides of the inequality above, and letting A go to 0, one 
obtains 

i-T 



E 



f e-^«(F(u(i))-Fo(i),w(i))di 
Jo 



> 0. 
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Since w is arbitrary, we conclude that Fo{t) — F{u{t)). Thus the existence of a 
strong solution of the stochastic Navier-Stokes equation (3.2) has been proved. □ 

4. Convergence as e — > 

Wc will now study the asymptotic limit of e ^ 0. Let us consider the family of 
perturbed systems(depending on the positive parameter e). 

2 1 

dtu' - uAu' + J2 <Au" + -(Divu^)u^ + Vp' = f + g{t)dw{t) (4.1) 

i=l 

in L'^{0,T;M-\0)), 
eatp^ + Divu'' = in L'^{0,T;L^{0)), (4.2) 

with the initial conditions 

u^(0) = uo in L2(C>) and p^{0)=po in ^^{0). (4.3) 

This is a method to overcome the computational difficulties connected with the 
constraint "Divu = 0". The equations (4.1)-(4.2) are easier to approximate than 
the original stochastic Navier-Stokes equation as the constraint "Divu = 0" has 

been replaced by the evolution equation (4.2). 

Here we will show how the solutions of the perturbed problems converge to the 
solutions of the incompressible stochastic Navier-Stokes equation as £ ^ 0. The 
idea of the proof is similar to the deterministic case presented in Temam [19] . 
Let 

f GL2(0,r;L2(O)), geL2(0,T;4(L2(O))), uoGh^{0), poEh'^iO). 

Then we can write the above mentioned perturbed systems as Ito Stochastic equa- 
tions in variational form 

' d{u'{t), v) + {-lyAu'it), v) dt + S(u^ u^ v)di 

+ {Vp%t),v}dt={f,v)dt + {git),^)dwit), (4.4) 
. {ep^{t) + Diyu%t),q) =0, 

in (0,T), with the initial conditions 

(u^(0),v) = (uo,v) and if{0),q) = ipo,q), (4.5) 

for any v in the space Hj(e») and any q in L^{n;L^{0,T;l?{O))). 

Proposition 4.1. As e ^ 0, the solutions {u'^,p^} of the equations (4.4)-(4.5) 
converge to the solution u of the incom,pressible stochastic Navier-Stokes equaiion. 

Proof. First we should point out that the solutions {u^,p^} of the equations (4.4)- 
(4.5) satisfy the monotonicity property in Lemma 2.4, the energy equality (3.7) 

and the a priori estimates (3.8) and (3.10). By virtue of these a priori estimates 
and using the Banach-Alaoglu theorem, along a subsequence the approximations 
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{u^,p^} have the following limits: 

— >u weakly star in LP(f];L°°(0,T;L2(C)))), (4.6) 
weakly in L2(0; L2(0, T; HJ(0))), 
y/Ip' — > X weakly in L^(f2; L^(0, T; L^(0))). (4.7) 
Let us denote 

F(u^) := Au' + B{u') - f := -uAu^ + [{u' ■ V) + (1/2) Divu^]u^ - f, 

and 

F(u) := Au + B{u) - f := -i^Au + (u • V)u - f . 
Let be a C°° scalar function on [0, T] with 4>{T) = 0. Multiplying the equation 
(4.4) by ^(t), integrating in t and taking mathematical expectation, we obtain 



+ E 



E 



Jo 







r 


-E 













+ E 











for all V in Hj(0). 



Now passing to the limit in (4.7) we have in the sense of distribution, 



E[V-ei%,,)]^E[i%,)]. 



(4.8) 



(4.9) 



Hence in the same sense 



0. 



Then passing to the limit in the following equation 

E[{ep^{t)+BiYu'{t),q)] =0 

we get 

E[{Divu,q)] =0, V9eL2(n;L2(0,T;L2(O))), 

which implies that Divu = 0, almost everywhere and almost surely. 
Hence u e LP(n; L~(0, T; I.^{0))) fl L^{n; L^{0, T; Hj(e>))). 
Now it is clear that 

(Vp^ v^(i)) = (p^ ^(i) Div v) = 0, almost surely. (4.10) 

Now using the same Minty-Browder monotonicity argument used in Proposition 
3.4 we can show that in limit 

F{u^) — > F{u) weakly in L2(n;L2(0,T;e-i(e)))). (4.11) 

Using the results (4.10) and (4.11), in limit we have from equation (4.8) 

-E 



[ {u{t)M'{mt\+E\ [ {F{u{t)),v<j>{t))dt 
Jo '-Jo 



=(uo,v)<^(0), 
for all V in Hj(0). 
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This proves that u is a solution of the incompressible stochastic Navier-Stokes 
equation. □ 

Remark 4.2. If one considers the multiplicative noise a{t, u) of the type considered 
in the hypotheses (^.1 — A.3) in Sritharan and Sundar [16], then under these 
conditions same a priori estimates (3.8)-(3.10) hold. Thus C7^(.,u'^) — » S weakly 
in L^(0; L^(0, T; Lq)), where Lq denote the space of linear operators S such that 
SQ^^"^ is a Hilbcrt- Schmidt operator from to and the norm on the space 
Lq is defined by \S\1^ = Ti{SQS*) and Q is a trace class operator. Hence with 
the help of Minty-Browder monotonicity arguments the existence and uniqueness 
proofs go through and we can also establish the limit to incompressible flow. 



Acknowledgment. We thank the reviewer for the helpful comments. 
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